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Abstract. The software Trim offers implementations of tropical implicitization and 
tropical elimination, as developed by Tevelev and the authors. Given a polynomial map 
with generic coefBcients, Trim computes the tropical variety of the image. When the 
image is a hypersurface, the output is the Newton polytope of the defining polynomial. 
Trim can thus be used to compute mixed fiber polytopes, including secondary polytopes. 
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1. Introduction. Implicitization is the problem of transforming a 
given parametric representation of an algebraic variety into its implicit 
representation as the zero set of polynomials. Most algorithms for elimi- 
nation and implicitization are based on multivariate resultants or Grobner 
bases, but current implementations of these methods arc often too slow. 
When the variety is a hypersurface, the coefhcients of the implicit equation 
can also be computed by way of numerical linear algebra |3l [7j , provided 
the Newton polytope of that implicit equation can be predicted a priori. 

The problem of predicting the Newton polytope was recently solved 
independently by three sets of authors, namely, by Emiris, Konaxis and 
Palios [5], Esterov and Khovanskii |13j . and in our joint papers with Tevelev 
[181 119| . A main conclusion of these papers can be summarized as follows: 
The Newton polytope of the implicit equation is a mixed fiber polytope. 

The first objective of the present article is to explain this conclusion 
and to present the software package Trim for computing such mixed fiber 
polytopes. The name of our program stands for Tropical Implicitization^ 
and it underlines our view that the prediction of Newton polytopes is best 
understood within the larger context of tropical algebraic geometry. The 
general theory of tropical elimination developed in |18j unifies earlier results 
on discriminants [3] and on generic polynomial maps whose images can have 
any codimension [TH] . The second objective of this article is to explain the 
main results of tropical elimination theory and their implementation in 
Trim. Numerous hands-on examples will illustrate the use of the software. 
At various places we give precise pointers to [8] and [13] , so as to highlight 
similarities and differences among the different approaches to the subject. 

Our presentation is organized as follows. In Section 2 we start out with 
a quick guide to Trim by showing some simple computations. In Section 
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3 we explain mixed fiber polytopes. That exposition is self-contained and 
may be of independent interest to combinatorialists. In Section 4 we discuss 
the computation of mixed fiber polytopes in the context of elimination 
theory, and in Section 5 we show how the tropical implicitization problem 
is solved in Trim. Theorem 5.1 expresses the Newton polytope of the 
implicit equation as a mixed fiber polytope. In Section 6 we present results 
in tropical geometry on which the development of Trim is based, and we 
explain various details concerning our algorithms and their implementation. 

2. How to use Trim. The first step is to download Trim from the 
following website which contains information for installation in Linux: 

jhttp : / /math .mit . edu/~jyu/TrIm, 

Trim is a collection of C-|— |- programs which are glued together and in- 
tegrated with the external software poljrmake [9] using perl scripts. The 
language perl was chosen for ease of interfacing between various programs. 

The fundamental problem in tropical implicitization is to compute the 
Newton polytope of a hypcrsurface which is parametrized by Laurent poly- 
nomials with sufficiently generic coefficients. As an example we consider 
the following three Laurent polynomials in two unknowns x and y with 
sufficiently generic coefficients ai , a2 , as , /3i , /32 , /^s , 7i , 72 , 73 : 

1 

u = ai ■ „ „ + a2 ■ X + as ■ xy 

9 1 

V = I3i-x^ + f32-y + Ps- - 

X 

2 1 1 

w = 71 • 2/ + 72 h 73 • -■ 

xy y 

We seek the unique (up to scaling) irreducible polynomial F{u^ v, w) which 
vanishes on the image of the corresponding morphism (C*)^ — > C'^. Using 
our software Trim, the Newton polytope of the polynomial F{u,v,w) can 
be computed as follows. We first create a file input with the contents 

[x,y] 

[x~(-2)*y"(-2) + X + x*y, 
x~2 + y + x~(-l) , 
y-2 + x-(-l)*y-(-l) + y-(-l)] 

Here the coefficients are suppressed: they are tacitly assumed to be generic. 
We next run a perl script using the command ./Trlm.prl input. The 
output produced by this program call is quite long. It includes the lines 

VERTICES 
1 9 2 
1 9 2 
1 9 
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1 9 
1 6 6 
1 2 2 8 
1 6 6 
1 2 8 2 
1 6 6 
1 
1 9 
1 2 9 
1 8 2 2 

Ignoring the initial 1, this hst consists of 13 lattice points in M."^, and these 
are precisely the vertices of the Newton polytope of F{u, v, w). The above 
output format is compatible with the polyhedral software Polymake [S]. 
We find that the Newton polytope has 10 facets, 21 edges, and 13 vertices. 
Further down in the output. Trim prints a list of all lattice points in the 
Newton polytope, and it ends by telling us the number of lattice points: 

N_LATTICE_PQINTS 
383 

Each of the 383 lattice points k) represents a monomial u^v^w^ which 
might occur with non-zero coefficient in the expansion of Fiu, v, w). Hence, 
to recover the coefhcients of F{u, v, w) we must solve a linear system of 382 
equations with 383 unknowns. Interestingly, in this example, 39 of the 383 
monomials always have coefficient zero in F{u, v, w). Even when ai, . . . , 73 
are completely generic, the number of monomials in F{u, v, w) is only 344. 

The command ./Trim.prl implements a certain algorithm, to be de- 
scribed in the next sections, whose input consists of n lattice polytopes 
in K"~^ and whose output consists of one lattice polytope in K". In our 
example, with n ~ 3, the input consists of three triangles and the output 
consisted of a three-dimensional polytope. These are depicted in Figure 1. 

The program also works in higher dimensions but the running time 
quickly increases. For instance, consider the hypersurface in represented 
by the following four Laurent polynomials in x, y, z, written in Trim format: 

[x,y,z] 

[x*y + z + 1, 

x*z + y + 1, 

y*z + X + 1, 

x~3 + y~5 + z"7] 

It takes Trim a few moments to inform us that the Newton polytope of 
this hypersurface has 40 vertices and contains precisely 5026 lattice points. 
The /-vector of this four-dimensional polytope equals (40, 111, 103, 32). 
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Output 



Fig. 1. Tropical implicitization constructs the three-dimensional Newton poly tope 
of a parametrized surface from the three Newton polygons of the given parametrization 

Remark 2.1. The examples above may serve as illustrations for the 
results in the papers [8] and [13]. Emiris, Konaxis and Palios [8] place the 
emphasis on computational complexity, they present a precise formula for 
plane parametric curves, and they allow for the map to given by rational 
functions. Esterov and Khovanskii develop a general theory of polyhedral 
elimination, which parallels the tropical approach in |18] . and which in- 
cludes implicitization as a very special case. A formula for the leading 
coefficients of the implicit equation is given in [8l §4]. This formula is cur- 
rently not implemented in Trim but it could be added in a future version. 

What distinguishes Trim from the approaches in [8] and [13] is the 
command TrCI which computes the tropical variety of a generic complete 
intersection. The relevant mathematics will be reviewed in Section 6. This 
command is one of the ingredients in the implementation of tropical implic- 
itization. The input again consists of m Laurent polynomials in n variables 
whose coefficients are tacitly assumed to be generic, or, equivalently, of m 
lattice polytopes in n-space. Here it is assumed that m < n. If equality 
holds then the program simply computes the mixed volume of the given 
polytopes. As an example, delete the last line from the previous input file: 

[x,y ,z] 

[x*y + z + 1, 

x*z + y + 1, 

y*z + X + 1] 



4 



The command . /TrCI . prl input computes the mixed volume of the three 
given lattice polytopes in R^. Here the given polytopes are triangles. The 
last line in the output shows that their mixed volume equals five. 

Now repeat the experiment with the input file input as follows: 
[x,y,z] 

[x*y + z + 1, x*z + y + 1] 

The output is a one-dimensional tropical variety given by five rays in M.^: 

DIM 
1 

RAYS 

-1 -1 

1 

10 

10 
-111 

MAXIMAL_CONES 



1 

2 

3 

4 

MULTIPLICITIES 

2 

1 

1 

1 

i 

Note that the first ray, here indexed by 0, has multiplicity two. This scaling 
ensures that the sum of the five RAYS equals the zero vector (0, 0, 0). 

For a more interesting example, let us tropicalize the complete inter- 
section of two generic hypersurfaccs in C^. We prepare input as follows: 

[a,b, c ,d, e] 

[a*b + b*c + c*d + d*e + a*e + 1, 
a*b*c + b+c*d + c*d*e + d*e*a + e*a*b] 

When applied to these two polynomials in five unknowns, the command 
./TrCI. prl input produces a three-dimensional fan in R^. This fan has 
26 rays and it has 60 maximal cones. Each maximal cone is the cone over 
a triangle or a quadrangle, and it has multiplicity one. The rays are 
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RAYS 

-110 1 
-111-13 
10 1 
-13-111 

The rays are labeled 0,1,..., 25, in the order in which they were printed. 
The maximal (three-dimensional) cones appear output in the format 

MAXIMAL_CONES 

12 3 

1 7 10 

1 12 

3 4 7 

3 12 

7 12 

It is instructive to compute the tropical intersection of two generic 
hypersurfaces with the same support. For example, consider the input file 

[x,y ,z] 

[l+x+y+z+ x*y + x*z + y*z + x*y*z, 
l+x+y+z+ x*y + x*z + y*z + x*y*z] 

As before, the reader should imagine that the coefficients are generic ra- 
tional numbers instead of one's. The tropical complete intersection de- 
termined by these two equations consists of the six rays normal to the six 
facets of the given three-dimensional cube. The same output would be pro- 
duced by Jensen's software GFan [11 [12], which computes arbitrary tropical 
varieties, provided we input the two equations with generic coefficients. 

3. Mixed fiber polytopes. We now describe the construction of 
mixed fiber polytopes. These generalize ordinary fiber polytopes [1], and 
hence they generalize secondary polytopes [TU]. The existence of mixed 
fiber polytopes was predicted by McDonald [T3| and Michiels and Cools (TB] 
in the context of polynomial systems solving. They were first constructed 
by McMuUen [TS] , and later independently by Esterov and Khovanskii [T3| . 

The presentation in this section is written entirely in the language of 
combinatorial geometry, and it should be of independent interest to some 
of the readers of Ziegler's text book [201. There are no polynomials or 
varieties in this section, neither classical nor tropical. The connection to 
elimination and tropical geometry will be explained in subsequent sections. 

Consider a linear map tt : — > and a p-dimensional polytope 
P C M.P whose image Q = 7r(P) is a q-dimensional polytope in R'^. li x is 
any point in the interior of Q then its fiber Tr~^{x) n P is a polytope of 
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dimension p — q. The fiber polytope is defined as the Minkowski integral 

S]^(P) = / {TT-\x)nP)dx. (3.1) 
Jq 

It was shown in [1] that this integral defines a polytope of dimension p — q. 
The fiber polytope (P) lies in an affine subspace of W which is a par- 
allel translate of kernel(7r). Billera and Sturmfels [T] used the notation 
E(P, Q) for the fiber polytope, and they showed that its faces are in bijec- 
tion with the coherent polyhedral subdivisions of Q which are induced from 
the boundary of P. We here prefer the notation S^(P) over the notation 
E(P, Q), so as to highlight the dependence on tt for fixed P and varying tt. 

Example 3.1. Let p = 3 and take P to be the standard 3-cube 

P = conv{(000),(001),(010),(011),(100),(101),(110),(lll)}. 

We also set q = 1 and we fix the linear map 

TT : R3 R\ {u,v,w) n- M + 2w + 3w. 

Then Q = tt{P) is the line segment [0,6]. For < a; < 6, each fiber 
Tr~^{x) n P is either a triangle, a quadrangle or a pentagon. Since the 
fibers have a fixed normal fan over each open segment (i, i + 1), we find 

s.(p) = ^ r^V'(-^)nP)dx = ^(^-i(, + l)np). 

Hence the fiber polygon is really just the Minkowski sum of two triangles, 
two quadrangles and two pentagon, and this turns out to be a hexagon: 

S^(P) = conv{(l, 10, 5), (1, 4, 9), (5, 2, 9), (11, 2, 7), (11, 8, 3), (7, 10, 3)}. 

In the next section we shall demonstrate how Trim can be used to compute 
fiber poly topes. The output produced will be the planar hexagon which is 
gotten from the coordinates above by applying the linear map (u,v,w) i— ^ 
— 3, w + ui — 9). Hence Trim produces the following coordinatization: 

I],(P) = conv{(2,6), (6,4), (6,2), (4,0), (0,2), (0,4) }. (3.2) 

It is no big news to polytope aficionados that the fiber polygon of the 3- 
cube is a hexagon. Indeed, by [501 Example 9.8], the fiber polytope obtained 
by projecting the p-dimensional cube onto a line is the permutohedron of 
dimension p — I. For p = 3 the vertices of the hexagon 'S^^^P) correspond 
to the six monotone edge paths on the 3-cube from (000) to (111). □ 

As a special case of the construction of fiber polytopes we get the 
secondary polytopes. Suppose that P is a polytope with n vertices in W 
and let A denote the standard (n — l)-simplex in R". There exists a linear 
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map p : R" — s- such that p(A) = P, and this hnear map is unique if 
we prescribe a bijection from the vertices of A onto the vertices of P. The 
polytope Ep(A) is called the secondary polytope of P; see |201 Definition 
9.9]. Secondary polytopes were first introduced in an algebraic context by 
Gel'fand, Kapranov and Zelevinsky [10]. For example, if we take P to be 
the 3-dimensional cube as above, then the simplex A is 7-dimensional, and 
the secondary polytope I]p(A) is a 4-dimensional polytope with 74 vertices. 
These vertices are in bijection with the 74 triangulations of the 3-cube. 

A detailed introduction to triangulations and a range of methods for 
computing secondary polytopes can be found in the forthcoming book [5]. 
We note that the computation of fiber polytopes can in principle be reduced 
to the computation of secondary polytopes, by means of the formula 

S,(p(A)) = p(S],op(A)). (3.3) 

Here tt o p is the composition of the following two linear maps of polytopes: 

A-^P^Q 

The formula (|3.3p appears in [U Lemma 2.3] and in [20l Exercise 9.6]. 
The algorithm of Emiris et al. [H §4] for computing Newton polytopes of 
specialized resultants is based on a variant of p.3|) . Neither our software 
Trim nor the Estcrov-Khovanskii construction [T^ uses the formula p.3p . 

We now come to the main point of this section, namely, the construc- 
tion of mixed fiber polytopes. This is primarily due to McMullen |15| , but 
was rediscovered in the context of elimination theory by Khovanskii and 
Esterov [13l §3]. Wc fix a linear map tt : M?* — ;> as above, but we 
now consider a collection of c polytopes Pi,...,Pc in R''. Wc consider 
the Minkowski sum Pa = AiPi + • • • + XcPc where A = (Ai, . . . , Ac) is a 
parameter vector of unspecified positive real numbers. We shall assume 
that Pa is of full dimension p, but we do allow its summands Pi to be 
lower-dimensional. The image of Pa under the map tt is the (/-dimensional 
polytope 

7r(PA) = Xi-tt{Pi) + --- + X,-tt{P,). 

The following result concerns the fiber polytope from Pa onto 7r(PA). 

Theorem 3.2 ([13 [IS])- The fiber polytope T.t,{Px) depends poly- 
nomially on the parameter vector A. This polynomial is homogeneous of 
degree q + 1. Moreover, there exist unique polytopes Mi-^i^...i^ such that 

S.(AiPi + ----f AcPc) = J2 ^^■■■K^-M,,^,...,^. (3.4) 

iiH hic=9+l 

To appreciate this theorem, it helps to begin with the case c = 1. That 
corresponds to scaling the polytopes P and Q above by the same factor A. 
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This results in the Minkowski integral (|3.ip being scaled by the factor A*'''^. 
More generally, the coefficients of the pure powers A^^^ in the expansion 
p.4p are precisely the fiber polytopes of the individual Pj, that is, 

Mo,...,0.q+1.0,....0 = ^Tl{Pj)- 

On the other extreme, we may consider ii = ^2 = ■ • ■ = ic = 1, which is the 
term of interest for elimination theory. Of course, if all ij's are equal to one 
then the number c of polytopes Pj is one more than the dimension q of the 
image of tt. We now assume that this holds, i.e., we assume that c = q+1. 
We define the mixed fiber polytope to be the coefficient of the monomial 
A1A2 • • • Ac in the formula (|3.4|) . The mixed fiber polytope is denoted 

^^{Pi,P2,...,Pc) := Mn...i. (3.5) 

The smallest non-trivial case arises when p = 3, c = 2 and <? = 1, where we 
are projecting two polytopes Pi and P2 in M.^. Their mixed fiber polytope 
with respect to a linear form tt : M'^ — !• is the coefficient of A1A2 in 

E,(AiPi + A2P2) = A?-S,(Pi) + AiA2-S,(Pi,P2) + A2.E,(P2). 

The following is [HI Example 4.10]. It will be revisited in Example 14.21 
Example 3.3. Consider the following two tetrahedra in three-space: 

Pi = conv{0, 3ei, 3e2, 3e3} and P2 = conv{0, — 2ei, — 2e2, —263}. 

Their Minkowski sum Pi + P2 has 12 vertices. 24 edges and 14 facets. If 
we take tt : M'^ — > M-'^ to be the linear form (u, v,w) 1-^ u — 2v + w then the 
fiber polytope S^(Pi -I- P2) = M20 + Mu + M02 is a polygon with ten 
vertices. Its summands AI20 = T,t^{Pi) and M02 = S^(P2) are quadrangles, 
while the mixed fiber polytope Mu ~ T,t^{Pi,P2) is a hexagon. □ 

We remark that fiber polytopes arc special instances of mixed fiber 
polytopes. Suppose that Pi = P2 = • ■ • = Pc are all equal to the same 
fixed polytope P in Rp. Then the fiber polytope SttCPa) in p.4p equals 

S(AiPi + • ■ • + AcPc) = (Ai + --- + A,)'=-S,(P). 

Hence the fiber polytope S7r(P) is the mixed fiber polytope SttCP, • ■ ■ , P) 
scaled by a factor of 1/cl. Similarly, any of the coefficients in the expansion 
(|3.4p can be expressed as mixed fiber polytopes. Up to scaling, wc have 

■^iii2 - ic = ^7r ( Pl, ■ ■ ■ , Pi , P2, ■ — 1 p2 _, ■ ■ ■ , Pc, ■ — , Pc ) ■ 

21 times 12 times ic times 

In the next section we shall explain how mixed fiber polytopes, and hence 
also fiber polytopes and secondary polytopes, can be computed using Trim. 
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4. Elimination. Let /i,/2,...,/c £ C[xf ^, x^"'^, . . . , ^] be Lau- 
rent polynomials whose Newton polytopes are Pi,P2,. ■ ■ ,Pc C W, and 
suppose that the coefficients of the fi are generic. This means that 

where the coefficients Ci^a are assumed to be sufficiently generic non-zero 
complex numbers. The corresponding variety 

X = {we (e)P : /i(h) = /2(u) = --- = /cH = 0} 

is a complete intersection of codimcnsion c in the algebraic torus (C*)^. 

We set 7' = p — c + 1 and we fix an integer matrix A = (oy- ) of format 
rxp where the rows of A are assumed to be linearly independent. We also 
let TT : RP R'^"^ be any linear map whose kernel equals the row space of 
A. The matrix A induces the following monomial map: 

a:{c*Y ^{c*Y-'+\ (xl,...,xp)K^(^x;^^...,^x;-). (4.1) 

J=l J=l 

Let Y be the closure in (C*)p^^+^ of the image a{X). Then F is a hypcr- 
surface, and wc are interested in its Newton polytopc. By this we mean 
the Newton polytope of the irreducible equation of that hypcrsurfacc. 

Theorem 4.1 (Khovanskii and Esterov [13]). The Newton polytope 
of Y is affinely isomorphic to the mixed fiber polytope S]^(Pi, . . . , Pc). 

A proof of this result using tropical geometry is given in [l^- The 
computation of the hypersurface Y from the defining equations /i , . . . , /c 
of X is a key problem of elimination theory. Theorem 14.11 offers a tropical 
solution to this problem. It predicts the Newton polytope of Y. This 
information is useful for symbolic-numeric software. Knowing the Newton 
polytopes reduces computing the equation of Y to linear algebra. 

The numerical mathematics of this linear algebra problem is interesting 
and challenging, as seen in [5] and confirmed by the experiments reported 
in [HJ §5.2]. We hope that our software Trim will eventually be integrated 
with software exact linear algebra or numerical linear algebra (e.g. LAPack). 
Such a combination would have the potential of becoming a useful tool for 
practitioners of non-linear computational geometry. 

In what follows, we demonstrate how Trim computes the Newton poly- 
tope of Y and hence the mixed fiber polytope . . . ,Pc)- The input 
consists of the polytopes Pi, . . . , Pc and the matrix A. The map tt is tacitly 
understood as the map from MP onto the cokernel of the transpose of A. 

Example 4.2. Let p = 3,c = 2 and consider [THl Example 1.3]. Here 
the variety X is the curve in (C*)'^ defined by the two Laurent polynomials 

/i = ai'x^ + a2X2 + ctsxl + and /2 = /3ix^'^ + (32X2'^ + PsX^'^ + fi^. 
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We seek to compute the Newton polygon of the image curve Y in (C*)^ 



[xl,x2,x3] 

[xl-3 + x2~3 + x3-3 + 1, xl~(-2)+x2~(-2)+x3~(-2)+l] 
We also prepare a second input file A. matrix as follows: 

LINEAR_MAP 
111 

12 

We now execute the following two commands in Trim: 

./TrCI.prl input > fan 

. /project .prl fcoi A. matrix 

The output we obtain is the Newton polygon of the curve Y: 

VERTICES 

1 36 
1 36 
1 30 12 
1 18 12 
1 6 24 
1 18 24 

This hexagon coincides with the hexagon in [181 Examples 1.3 and 4.10]. 
It is isomorphic to the mixed fiber polytope I]^(Pi,P2) hr Example 13.31 □ 

We may use Trim to compute arbitrary fiber polytopes. For example, 
to carry out the computation of Example 13.11 we prepare input as 

[x , y , z] 

[l+x+y+z+ x*y + x*z + y*z + x*y*z, 
l+x+y+z+ x*y + x*z + y*z + x*y*z] 

and A. matrix as 

LINEAR.MAP 

11-1 

2-10 

The two commands above now produce the hexagon in p.2p . Our next 
example shows how to compute secondary polytopes using Trim. 



Example 4.3. Following [20l Example 9.11], we consider the hexagon 
with vertices (0, 0), (1, 1), (2, 4), (3, 9), (4, 16), (5, 25). This hexagon is rep- 
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resented in Trim by the following file A.matrix. The rows of this matrix 
span the linear relations among the five non-zero vertices of the hexagon: 

LINEAR_MAP 

3 -3100 

8 -6010 
15 -10 1 

On the file input we take three copies of the standard 5-simplex: 
[a,b,c,d,e] 

[a+b+c+d+e+1 , a+b+c+d+e+1, a+b+c+d+e+1] 

Running our two commands, we obtain a 3-dimensional polytope with 14 
vertices, 21 edges and 9 facets. That polytope is the associahedron [20]. □ 

We close this section with another application of tropical elimination. 

Example 4.4. For two subvarieties Xi and X2 of (C*)" we define their 
coordinate-wise product Xi ★ X2 to be the closure of the set of all points 
(wiwi,...,u„t;„) where (ui,...,m„) e Xi and («!,...,«„) £ X2. The 
expected dimension of Xi -k X2 is the sum of the dimensions of Xi and X2, 
so we can expect Xi -k X2 to be a hypersurface when dim(Xi) + dim(X2) = 
71 — 1. Assuming that Xi and X2 are generic complete intersections then 
the Newton polytope of that hypersurface can be computed using Trim as 
follows. Let p = 2n and define X as the direct product Xi x X2. Then 
Xi ★ X2 is the image of X under the monomial map 

a : (C*)^" (C*)", (Ui, ...,Un,Vi,...,Vn)'-^ (uivi,. . . ,u„u„). 

Here is an example where Xi and X2 are curves in three-dimensional space 
(n = 3). The two input curves are specified on the file input as follows: 

[ul,u2,u3, vl,v2,v3] 
[ul + u2 + u3 + 1, 

ul*u2 + ul*u3 + u2*u3 + ul + u2 + u3, 
vl*v2 + vl*v3 + v2*v3 + vl + v2 + v3 + 1, 
vl*v2*v3 + vl*v2 + vl*v3 + v2*v3 + vl + v2 + v3] 

The multiplication map a : (C*)^x(C*)^ (C*)^ is specified on A.matrix: 

LINEAR_MAP 
10 10 
10 10 
1 1 

The image of Xi x X2 under the map a is the surface Xi * X2 ■ We find 
that the Newton polytope of this surface has ten vertices and seven facets: 
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VERTICES 
18 4 
10 8 4 
10 8 
10 8 
14 8 
14 8 
10 4 8 
18 
10 
18 4 

FACETS 

128 -16 

32 
128 -16 

32 

64 
128 -16 
192 -16 -16 -16 

5. Implicitization. Implicitization is a special case of elimination. 
Suppose we are given n Laurent polynomials 171, . . . , g„ in C[tf^, . . . , t^^-^] 
which have Newton polytopes Qi, . . . , Qn C M"~^ and whose coefficients 
are generic complex numbers. These data defines the morphism 

g : iC*r-' ^ {C*y\ t ^ (giW,...,g„_i(0). (5.1) 

Under mild hypotheses, the closure of the image of g is a hypersurface 
Y in (C*)". Our problem is to compute the Newton polytope of this 
hypersurface. A first example of how this is done in Trim was shown in the 
beginning of Section 2, and more examples will be featured in this section. 

The problem of implicitization is reduced to the elimination compu- 
tation in the previous section as follows. We introduce n new variables 
2/1, . . . , ?/„ and we consider the following n auxiliary Laurent polynomials: 

hix) = 9i{t) - yi , h{x) = g-iit) - 2/2 , • • ■ , fn{x) = gn{t) - Vn- (5.2) 

Here we set p = 2n — 1 and (xi, . . . , Xp) = (ii, . . . , t„-i, yi, . . . , y„) so as to 
match the earlier notation. The subvariety of (C*)^ = (C*)"-i x (C*)" 
defined by /i, ...,/„ is a generic complete intersection of codimcnsion n, 
namely, it is the graph of the map g. The image of g is obtained by pro- 
jecting the variety {/i = •■• = /« = 0} onto the last n coordinates. This 
projection is the monomial map a specified by the n x p-matrix A = (0 /) 
where is the n x (n— 1) matrix of zeroes and / is the n x n identity matrix. 

This shows that we can solve the implicitization problem by doing the 
same calculation as in the previous section. Since that calculation is a main 
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application of Trim, we have hard- wired it in the command . /Trim.prl. 
Here is an example that illustrates the advantage of using tropical implic- 
itization in analyzing parametric surfaces of high degree in three-space. 

Example 5.1. Consider the parametric surface specified by the input 
[x,y] 

[x~7*y"2 + x*y + x~2*y~7 + 1, 
x"8*y"8 + x~3*y~4 + x~4*y~3 + 1, 
x"6*y + x*y~6 + x~3*y~2 + x~2*y~3 + x + y] 

Using the technique shown in Section 2, we learn in a few seconds to learn 
that the irreducible equation of this surface has degree 90. The command 
. /Trim.prl input reveals that its Newton polytope has six vertices 

VERTICES 
1 80 
1 45 
1 80 
1 10 80 
1 
1 28 54 

This polytope also has six facets, namely four triangles and two quadran- 
gles. The expected number of monomials in the implicit equation equals 

N_LATTICE_POINTS 
62778 

At this point the user can make an informed choice as to whether she wishes 
to attempt solving for the coefficients using numerical linear algebra. □ 

Returning to our polyhedral discussion in Section 3, we next give a 
conceptual formula for the Newton polytope of the implicit equation as 
a mixed fiber polytope. The given input is a list of n lattice polytopes 
Qi, Q2, . . . , Qn in M"^^. Taking the direct product of with the space 
M" with standard basis {ei, 62, . . . , e„}, we consider the auxiliary polytopes 

conv(QiU{ei}), conv(Q2 U {ez}) , conv(Q„ U {e„}) C M''^^ x M", 

where conv denotes the convex hull. These are the Newton polytopes of 
the equations /i, /2, • . • , /« in (|5.2p . We now define tt to be the projection 
onto the first n — 1 coordinates 

TT : R""i X M" R""i , (ui, . . . ,Un-i,Vi,V2, ■■.,vn) 1-^ (mi, . . . ,U„-l). 

The following result is an immediate corollary to Theorem 14. II We propose 
that it be named the Fundamental Theorem of Tropical Implicitization. 
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Theorem 5.2. The Newton polytope of an irreducible hypersurface in 
(C*)" which is parametrically represented by generic Laurent polynomials 
with given Newton polytopes Qi, . . . , Qn equals the mixed fiber polytope 

(conv((3i U {ei}), conv((32 U {62}), . . . , conv((5„ U {e„})) . (5.3) 

This theorem is the geometric characterization of the Newton poly- 
tope of the imphcit equation, and it summarizes the essence of the recent 
progress obtained by Emiris, Konaxis and Pahos [8], Esterov and Khovan- 
skii [131, and Sturmfels, Tevelev and Yu [IHlIini- Our implementation of in 
Trim computes the mixed fiber polytope (|5.3p for any given Qi,Q2, ■ ■ ■, Qn, 
and it suggests that the Fundamental Theorem of Tropical Implicitization 
will be a tool of considerable practical value for computational algebra. 

Example 5.3. We consider a threefold in which is parametrically 
represented by four trivariate polynomials. On the file input we write 

[x,y ,z] 

[x + y + z + 1, 
x"2*z + y~2*x + z~2*y + 1, 
x"2*y + y~2*z + z~2*x + 1, 
x*y + x*z + y*z + x + y + z] 

The Newton polytope of this threefold has the f-vector (8, 16, 14, 6), and it 
contains precisely 619 lattice points. The eight vertices among them are 

VERTICES 
1 15 
1 6 
1 9 
1 6 
1 
1 12 3 
1 9 3 
1 9 3 

This four-dimensional polytope is the mixed fiber polytope (|5.3p for the 
tetrahedra Qi, (52, and the octahedron Q/^ specified in the file input. □ 

In (|5.ip we assumed that the gi{t) are Laurent polynomials but this 
hypothesis can be relaxed to other settings discussed in [HI [13] . In partic- 
ular, the Fundamental Theorem of Tropical Implicitization extends to the 
case when the gi{t) are rational functions. Here is how this works in Trim. 

Example 5.4. Let ai, . . . , as and . . . , /^s be general complex num- 
bers and consider the plane curve which has the rational paramctrization 



ait^ + a2t + as 



and y 



(5.4) 



X = 



ait'^ -t- as 
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This curve appears in [H Example 4.7]. The input for Trim is as follows: 
[ t, X, y ] 

[ t~3 + t + 1 + x*t"2 + X , 
t~4 + t~3 + 1 + y*t~2 + y ] 

The equation of the plane curve is gotten by eliminating the unknown 
t from the two equations (|5.4p . Tropical elimination using Trim predicts 
that the Newton polygon of that plane curve is the following pentagon: 

POINTS 
14 2 
10 3 
12 3 
10 
14 

This prediction is the correct Newton polygon for generic coefficients ai 
and In particular, generically the pairs of coefficients (a4,a5) and 
{Pii P^) in the denominators arc distinct. If we assume, however, that 
the denominators arc the same, then the correct Newton polytope is a 
quadrangle inside this pentagon, as computed in [5J Example 5.6]. □ 

6. Tropical varieties. We now explain the mathematics on which 
Trim is based. The key idea is to embed the study of Newton polytopes 
into the context of tropical geometry [2j |4j [HI [19] . Let / be any ideal in 
the Laurent polynomial ring iC[xf^, . . . ,x^^]. Then its tropical variety is 

T{I) = {w E W : in^(/) does not contain a monomial}. 

Here inu,(/) is the ideal of C[a;f^, . . . ,2;^^] which is generated by the w- 
initial forms of all elements in /. The set T{I) can be given the structure of 
a polyhedral fan, for instance, by restricting the Grobncr fan of any homog- 
enization of /. A point w in T{I) is called regular if it lies in the interior 
of a maximal cone in some fan structure on T{I). Every regular point w 
naturally comes with a multiplicity m.^,, which is a positive integer. We 
can define the sum of multiplicities of all minimal associate primes 

of the initial ideal in^ (/) . The multiplicities TOu, on 7~(/) are independent 
of the fan structure and they satisfy the balancing condition |18[ Def. 3.3]. 

If / is a principal ideal, generated by one Laurent polynomial f{x), 
then T{I) is the union of all codimension one cones in the normal fan of 
the Newton polytope P oi f{x). A point w € T{I) is regular if and only 
if w supports an edge of P, and is the lattice length of that edge. It 
is important to note that the polytope P can be reconstructed uniquely, 
up to translation, from the tropical liypersurfacc T{I) together with its 
multiplicities m^. The following Trim example shows how to go back and 
forth between the Newton polytope P and its tropical hypersurface 7~(/). 

16 



Example 6.1. We write the following polynomial onto the file poly: 



[ X, y, z ] 

[ X + y + z + x~2*y~2 + x~2*z~2 + y"2*z~2 ] 

The command ./TrCl.prl poly > fan writes the tropical surface defined 
by this polynomial onto a file fan. That output file starts out like this: 

AMBIENT_DIM 
3 

DIM 
2 



The tropical surface consists of 12 two-dimensional cones on 8 rays in R^. 
Three of the 12 cones have multiplicity two, while the others have multi- 
plicity one. Combinatorially, this surface is the edge graph of the 3-cube. 
The Newton polytope P is an octahedron, and it can be recovered from the 
data on fan after we place the 3x3 identity matrix in the file A. matrix: 

LINEAR_MAP 
10 
10 
1 

Our familiar command 
./project .prl fan A. matrix 

now reproduces the Newton octahedron P in the familiar format: 

VERTICES 
12 2 
10 2 2 
10 10 
12 2 
110 
10 1 

Note that the edge lengths of P are the multiplicities on the tropical surface. 
□ 

The implementation of the command project. prl is based on the 
formula given in |19l Theorem 5.2]. See [H §2] for a more general version. 
This result translates into an algorithm for Trim which can be described as 
follows. Given a generic vector w G K'^ such that face^(P) is a vertex v, the 
i*^ coordinate Wj is the number of intersections, counted with multiplicities, 
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of the ray u; + M>oei with the tropical variety. Here, the muhiphcity of the 
intersection with a cone F is the multiphcity of F times the absolute value 
of the I*'' coordinate of the primitive normal vector to the cone F. 

Intuitively, what we are doing in our software is the following. We 
wish to determine the coordinates of the extreme vertex v = faceu,(P) of 
a polytope P in a given direction w. Our polytope is placed so that it 
lies in the positive orthant and touches all the coordinate hyperplanes. To 
compute the i"" coordinate of the vertex v, we can walk from v toward the 
i*'' hyperplane along the edges of P, while keeping track of the edge lengths 
in the z*'' direction. A systematic way to carry out the walk is to follow the 
edges whose inner normal cone intersect the ray u; + M>oei. Recall that the 
multiplicity of a codimension one normal cone is the lattice length of the 
corresponding edge. Using this subroutine for computing extreme vertices, 
the whole polytope is now constructed using the method of Huggins [TT] . 

To compute the tropical variety ?'(/) for an arbitrary ideal I one can 
use the Grobner-based software GFan due to Jensen [21 [12] . Our polyhedral 
software Trim performs the same computation faster when the generators 
of / are Laurent polynomials /i , /2 , ■ ■ • , /c that are generic relative to their 
Newton polytopcs Pi, P2, . . . , Pc. It implements the following combinato- 
rial formula for the tropical variety T{I) of the complete intersection /. 

Theorem 6.2. The tropical variety T{I) is supported on a subfan of 
the normal fan of the Minkowski sum X]i=i ^i- ^ point w is in T{I) if and 
only if the polytope facCn,(^^g j Pj) has dimension > \J\ for J C {1, . . . , c}. 
The multiplicity of T{I) at a regular point w is the mixed volume 

m.u] = mixed volume(faceu,(Pi), face^(P2), face^(Pc)) , (6.1) 

where we normalize volume respect to the affine lattice parallel to X^jeJ ^i' 

For a proof of this theorem see [IHl §4] . We already saw some examples 
in the second half of Section 2. Here is one more such illustration: 

Example 6.3. We consider the generic complete intersection of codi- 
mension three in six-dimensional space (C*)"^ given by the polynomials 

[a,b,c,d,e,f] 

[a*b*c*d*e*f +a+b+c+d+e+f+l, 

a*b + b*c + c*d + d*e + e*f + f*a, 

a + b + c + d + e + f + 1] 

The application of . /TrCI . prl to this input internally constructs the cor- 
responding six-dimensional polytope Pi -I- P2 -I- P3. It lists all facets and 
their normal vectors, it generates all three-dimensional cones in the nor- 
mal fan, and it picks a representative vector w in the relative interior of 
each such cone. The three polytopes faceu,(Pi), faceuj(P2) and faceu,(P3) 
are translated to lie in the same three-dimensional space, and their mixed 
volume is computed. If is positive then Trim outputs the rays of 
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that cone and the mixed volume . The output reveals that this tropical 
threefold in consists of 117 three-dimensional cones on 22 rays. □ 

In our implementation of Theorem 16.21 in Trim, the Minkowski sum 
is computed using the iB4e library |llj . enumerating the d-dimensional 
cones is done by Polymake [^j, the mixed volumes are computed using the 
mixed volume library [6], and integer linear algebra for lattice indices 
is done using NTL [17]. What oversees all the computations is the perl 
script TrCl.prl. The output format is consistent with that of the current 
version of Gf an [12j and is intended to interface with the polyhedral software 
Polymake [3] when it supports polyhedral complexes and fans in the future. 

Elimination theory is concerned with computing the image of an alge- 
braic variety whose ideal / we know under a morphism a : (C*)^ — ?> (C*)''. 
We write j3 for the corresponding homomorphism from the Laurent poly- 
nomial ring in r unknowns to the Laurent polynomial ring in p unknowns. 
Then J = /3~^{I) is the ideal of the image variety, and, ideally, we would 
like to find generators for J. That problem is too hard, and what we do 
instead is to apply tropical elimination theory as follows. We assume that 
a is a monomial map, specified by an r x p integer matrix A = (a^ ) as 
in (|4.1|) . Rather than computing the variety of J from the variety of /, we 
instead compute the tropical variety 7~(J) from the tropical variety T(/). 
This is done by the following theorem which characterizes the multiplicities. 

Theorem 6.4. The tropical variety T{J) equals the image of 
under the linear map A. // the monomial map a induces a generically 
finite morphism of degree S from the variety of I onto the variety of J then 
the multiplicity of T{J) at a regular point w is computed by the formula 

my, = i • y mu • index(L„ nZ"^ : A(L,„ nZ^)). (6.2) 

V 

The sum is over all points v in T{I) with A.v — w. We assume that the 
number of these points is finite, they are all regular in T{I), and L„ is the 
linear span of a neighborhood of v in T{I), and similarly for w G T{J)- 

The formula (j6.2p can be regarded as a push-forward formula in in- 
tersection theory on toric varieties, and it constitutes the workhorse inside 
the Trim command project.prl. When the tropical variety T(J) has 
codimension one, then that tropical hypcrsurface determines the Newton 
polytope of the generator of J . The transformation from tropical hyper- 
surface to mixed fiber polytope was behind all our earlier examples. That 
transformation was shown explicitly for an octahedron in Example 16.11 

In all our examples so far, the ideal / was tacitly assumed to be a 
generic complete intersection, and Theorem 16. 21 was used to determine the 
tropical variety ?"(/) and the multiplicities m^. In other words, the com- 
mand TrCl furnished the ingredients for the formula (|6.2p . In particular, 
then / is the ideal of the graph of a morphism, as in Section 5, then Theo- 
rem [^13] specializes to the formula for tropical implicitization given in |19| . 
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It is important to note, however, that Theorem 16 .41 apphes to any ideal 
I whose tropical variety happens to be known, even if / is not a generic 
complete intersection. For instance, 7"(/) might be the output of a GFan 
computation, or it might be one of the special tropical varieties which have 
already been described in the literature. Any tropical variety with known 
multiplicities can serve as the input to the Trim command project.prl. 
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